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A universal description of ultraslow glass dynamics 
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The dynamics of glass is of importance in materials science but its nature has not yet been 
fully understood. Here we report that a verification of the temperature dependencies of the 
primary relaxation time or viscosity in the ultraslowing/ultraviscous domain of glass-forming 
systems can be carried out via the analysis of the inverse of the Dyre-Olsen temperature 
index. The subsequent analysis of experimental data indicates the possibility of the self- 
consistent description of glass-forming low-molecular-weight liquids, polymers, liquid crys- 
tals, orientationally disordered crystals and Ising spin-glass-like systems, as well as the 
prevalence of equations associated with the 'finite temperature divergence'. All these lead to 
a new formula for the configurational entropy in glass-forming systems. Furthermore, a link to 
the dominated local symmetry for a given glass former is identified here. Results obtained 
show a new relationship between the glass transition and critical phenomena. 
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The glass transition constitutes one of the outstanding 
problems of condensed matter physics 1 , with fundamental 
importance for many areas ranging from geophysics to 
modern food science and material engineering 2-5 . Over the last 
decades, an important progress in this field has been achieved, 
primarily because of the progress in the state-of-the-art 
experimental, theoretical and numerical insights. Hence, a 
conceptual breakthrough seems to emerge 6-18 . 

Recently, however, serious doubts regarding the most funda- 
mental experimental artifact, namely the parameterization of 
viscosity rj(T) or primary relaxation time t(T) on approaching the 
glass temperature (T g ), appeared 15 ' 19-21 . The importance of this 
issue is related to the fact that it is often associated with a 
hypothetical link between the glass transition and critical 
phenomena physics, that is, the thrilling possibility of a new 
great unification within condensed matter physics 14 ' 1 . Moreover, 
relations portraying the evolution of >j(T) or t(T) are often 
considered as a checkpoint for theoretical models 3 ' 13 . Decades of 
studies have led to the prevailing conviction that the ultimate 
parameterization is possible via the Vogel-Fulcher-Tammann 
(VFT) equation 3 ' 13 ' 22-24 : 



t(T) = t 0 exp 



DtTq 
T-T 0 



and rj(T) = r\ a exp 



DtTq 
T-T 0 



(1) 



where D T denotes the fragility strength coefficient and T 0 is the 
VFT singular temperature (T 0 <C T g and T> T g ). 

However, by following the compelling analysis of 43 sets of 
t(T) data for low-molecular- weight ultraviscous liquids, Hecksher 
et al. 15 recently concluded: 'There is no compelling evidence for 
the VFT prediction that the relaxation time diverges at a finite 
temperature [ . . . ] with a dynamic divergence of the VFT form 
lack a direct experimental basis .... we suggest that in the search 
for the correct theory for ultraviscous liquid dynamics, theories 
not predicting a dynamic divergence of the VFT form should be 
focused on'. McKenna 25 extended this outlook to the case of 
polymer glass formers. These findings were strengthened by the 
derivation by Mauro et al} 9 of the equation introduced 
empirically by Waterton 26 : 



t(T) = t 0 exp 



K 



r exp vr 



c 



The subsequent analysis of experimental data confirmed the 
advantage of the Waterton-Mauro (WM, called also MYEGA) 
equation without a 'finite temperature divergence' over the VFT 
one ii. 19-21,27 fh,^ me question regarding t(T) or r\(T) portrayal 
seemed to be clarified. However, the problem becomes more 
complex if glass-forming liquid crystals (LCs) and orientationally 
disordered crystals (ODICs; plastic crystals) are also considered. 
For these materials, the prevalence of the 'critical-like' equation 
with the clear 'finite temperature divergence' was clearly 
shown 28 ' 29 : 



(3) 



where T > T g and T c < T g : most often T g — T 0 « 30-40 K and 
T g - T c kW-15K. 

Surprisingly, the latter description seems to be applicable also 
for selected low-molecular-weight liquids and polymers 28 ' 29 . 
Worth mentioning is that such dynamics is also commonly 
accepted for spin-glass-like (SGL) systems 5 ' 28 . Remarkable is 
furthermore the approximate coincidence of power exponents in 
the mentioned systems, namely </>«9 for LC 28 ' 29 , </> = 9—15 in 
ODIC 28 ' 29 and 0 = 9-12 in SGL 3 ' 13 ' 28 . 

In this respect, it is important to point out the link between 
anomalous dynamics in supercooled glass-forming systems and 



the cooperatively rearranged regions, introduced by the Adam 
and Gibbs (AG) model 3 ' 3 ' 30 , leading to the following output 
relation 30 : 



t(T) = t 0 exp I ' — -I 



( AA/A 
\TSc) 



(4) 



where A/i defines that the free-energy barrier between 
cooperatively rearranged region and S c is the configurational 
entropy related to the difference between the entropy of the 
metastable disordered supercooled system and the stable crystal 
state. 

Angell et aZ. 31-33 noted that the AG equation can be converted 
into the VFT one by assuming: 



Sc(T)=S 0 



r-r c 

T 



(5) 



for T 0 = T K , where the latter denotes the Kauzmann temperature 34 . 

Although equation (5) became popular in the semi-phenom- 
enological analysis of supercooled glass formers 3 ' , some 
puzzling issues remain. Mauro 11 ' 19 has shown that the 
'divergence-less' equation (2) can also be derived from the AG 
dependence by omitting equation (5) for the configurational 
entropy. Moreover, the analysis of experimental data indicated 
that the coincidence between T 0 and T K may be limited 35 . 
Eckmann and Procaccia 36 showed that for a soft potential based 
model, the configurational entropy can stay positive for T > 0, 
concluding that the Kauzmann temperature may not exist. 

Particularly noteworthy is that the parameterization of t(T) or 
rj(T) via 'finite temperature divergence' relations was used as an 
important argument to support the hypothesis, which correlates 
the glass transition with a hidden phase transition below T g . For 
instance, Colby 16 suggested that equation (3), with exponent 
4> = 9, may be 'universal' for glass-forming low-molecular-weight 
liquids and polymers, although for some of them, a multiplicative 
factor exp(E/kT) should be included. The 'critical-like' divergence 
was associated with a hypothetical phase transition at T c < T g , 
linked to the correlation length £(T) = ? 0 (T- T c ) " v with 
v = 3/2 (refs 16,17). However, a subsequent thorough analysis 
of t(T) experimental data did not confirm the suggested superior- 
fitting features 28 ' 29 ' 37 . Notwithstanding, there are several low- 
molecular-weight and polymeric glass formers for which the 
behaviour suggested by Colby holds surprisingly well 28 ' 29 . 

It is worth mentioning that in 2010, Tanaka et aZ. 18 carried out 
an in-depth simulation analysis of six glass-forming systems and 
stated: 'Our results suggest a far more direct link than thought 
before between glass transition and critical phenomena. Indeed, 
the glass transition may be a new type of critical phenomenon 
where a structural order parameter is directly linked to slowness'. 
Tanaka et aZ. 18 linked the 'pretransitional' behaviour of the 
relaxation time to the Ising-like critical power-law divergence of 
the correlation length towards the ideal glass transition with the 
exponent v « 2/d, where d defines the spatial dimensionality. 
This led to the r VFT-type equation t(T) = x x exp[D r T 0 / 



(T-r 0 ) 



Tanaka et al. , further 



.explDHctTVcor 2 
suggested that glass transition may t>e encountered in the 'Ising 
universality class for systems with a non-conserved order 
parameter'. 

At this point, some basic concepts of critical phenomena 
physics are worth recalling 38 ' 39 . The pretransitional behaviour for 
different physical magnitudes {X(T}) is described by the relation 
X(T)cc\T — T c \^, where the values of critical exponents (1) are 
universal and depend solely on the dimensionality of the order 
parameter, the space dimensionality and the range of interactions. 
The symmetry is playing an essential role, whose changes are 
quantified by the order parameter. Microscopically, different 
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systems can be ascribed in universality classes. It is also worth to 
bear in mind that the physics of critical phenomena is considered 
one of the greatest 'unification' successes in condensed matter 
physics of the 20th century 2 ' 38 ' 39 . To the best of the authors' 
knowledge, no clear manifestation of the above fundamental 
features for the glass-transition phenomenon has been reported so 
far. This paper aims to fill this gap and resolve the Gordian knot of 
the t(T) or ii(T) evolution on approaching the glass transition. 

Results 

Fragility and the Dyre-Olsen temperature index. The Super- 
Arrhenius (SA) evolution of relaxation time or viscosity in the 
ultraslowing/ultraviscous domain can be generally expressed as 
z(T) = x oa exp[^(t)/RT], where £ a (T) = 7/ a (T) / R stands for 
the normalized apparent activation energy (AAE) and R denotes 
the gas constant ' . The departure from the simple Arrhenius 
behaviour with £ a = const is characterized by the fragility 
coefficient (m), introduced by Angell et al. 40,41 as the metric at 
the normalized plots log 10 T or log 10 >7 versus T g /T: 

_ \ rdlog,„ t(T) 
m = m v (T=T B ) = 



weight glass-forming ultraviscous liquids. To quantify how 



■mAT=T„) = 



. d(r g /r) . 

dl °gio'?( T ) 



d(r g /r) 



(6) 



Although the fragility coefficient m provides a metric for glass- 
forming liquids within a wide diversity of materials, it is not the 
only parameter of this type 3 ' 13 . A decade ago, Dyre-Olsen (DO) 42 
noted the formal importance of the Griineisen parameter (y) 43 : 



V d(Oi d In C0j 
(OjlTV ~ ~ dlnV 



(7) 



It indicates that if the restoring force acting on an atom is non- 
linear in displacement, the individual phonon frequency co, 
changes with the volume of crystal V, and parameter y reaches a 
transition from zero (linear displacement) to y > 0 (ref. 43). In the 
SA domain, the AAE increases with a decreasing temperature, 
allowing a dramatic increase from constant values to 
dE a (T)/dT > 0, hence giving rise to a non-linear behaviour in 
the log 10 T versus 1/T representation. This formal similarity 
allowed to propose the metric of the SA behaviour by the 
'Gruneisen-style' activation energy temperature index 42 : 

dlnEJT) , . 

The DO temperature index can serve as the alternative 
metric of fragility, which can be also presented as 
m=m(T=r g ) = (H;(r g )/T g )log loe , where ^(T) = HjR 
and HJT) = dlnT/d(l/T) denote the apparent activation 
enthalpy 4 

Assuming the hypothetical 'universal' value for the 
prefactor t 0 =10 s (refs 3,13,15), one obtains: 
m = log 10 [r(r g )/T 0 ] (lD 0 (r g _) + 1) = 16[/Do(r g ) + 1], for 
T(T g ) = 10 s and t 0 =10 s. Subsequently, taking into 
account the recent suggestion that 16<m<175 (ref. 45), one 
can estimate: 0 < 7 DO (T g ) < 10. It is noteworthy that in 1997, 
Simon et al. 46 indicated the upper fragility limit m = 214, yielding 
7 DO (T g ) = 12.375. Very recently, McKenna 47 indicated that if T 0 
is an ideal glass transition, the value of m evaluated at T= T 0 = T g 
would tend towards infinity for all materials. The same can be 
expected for 7 DO (T„). 

Heksher et al. showed that the analysis of I DO Cr) may 
represent a decisive tool for testing the validity of the given 
equation for fitting t(T) experimental data in low-molecular- 



exactly the VFT and the popular Bassler-Avramov' 



48,49 



(BA) 



t(T) = t 0 exp(j4/T D ) equations describe the data, they compared 
experimental and 'theoretical' values of 7 D o(T), namely 15 : 



T-Tn 



/BA( T ) =D _ 1 = (;onst 



(9) 



Derivations of Ioo(T) for equations (2) and (3) are presented in 
the Methods section. In the ultraslowing domain, near T g , they 
are as follows: 



^DO 



(T): 



C 
T 



T-T c 



(10) 



A glimpse on equations (9) and (10) shows that reciprocals of 
DO indexes yield a linear dependence I^q(T) = b + al, namely: 



ko(T) 



f(*) r 

i 

D- 1 ' 



- 1 -> (VFT) 

(BA) 
* (WM) 



(11) 



(Grit.; 



assuming 



14 s as an 'universal averaged' value of the prefactor. 
in the systems analysed in this paper 



The discussed evolutions of DO indexes are shown graphically 
in Fig. 1. Two characteristic features are worth mentioning, 
namely: 7jJq (T = 0) = b and 7jJ 0 ' (7c, T 0 ) = 0. It is worth stressing 
that while the physics associated with the fragility index m has 
been very broadly discussed 3 ' 13 ' 15,20 ' 31 ' 40-45 ' ' 51 , the meaning 
and properties of the DO index is poorly explored so far 15 ' 42 . 

For the estimation of 7 DO (7*), the evolution of EJJ) has to be 
determined. This, however, constitutes a notable issue. Hecksher 
et al} 5 calculated 7s a (T) from the SA relation 
To=10"" 
However. 

t 0 = 10" n s^-10" 16 s (refs 28,29,44). In the Methods section, 
we present a novel way of determining 74(7), avoiding the biasing 
influence of poorly estimated values of the prefactor t 0 for the 
given system. 

Experimental data analysis. This paper focuses on the ultraslow/ 
ultraviscous domain as in Hecksher et al} 5 , for x > t(T b ) « 
10 " 7 s or }] > )?(T B ) > 10 3 Poise, which usually occurs for 
T B — T g < 80K (refs 50-53). By heating above T B , the so-called 
dynamic cross-over occurs. The discussion on recent advances 
regarding this phenomenon can be found in Martinez-Garcia 
et al, 52 The analysis of experimental data presented below does 
not only cover the most often discussed supercooled low- 
molecular-weight liquids (L) and multimeric and polymeric (P) 
glass formers as in Hecksher et al} 5 , but also LCs, plastic crystals 
(ODICs) and SGL systems. The analysis is based on 16 selected 
sets of t(T) experimental data representing these types of glass- 
forming systems. However, to clarify the discussion, the key 
results are presented for a smaller number of characteristic 
systems. They are L including diethyl phthalate (DEP), glycerol, 
sorbitol, ethanol 28 ' 54 , P including EPON 828 and polystyrene 
(M w = 700) (refs 28,55), ODIC (plastic crystal) including C7-OH, 
C8-OH (ref. 53), CNc6 (ref. 57), CNadm 58 - 59 and NPaNPG 03 
(ref. 56), and SGL including ferrofluid liquid, with 5% solution of 
single-domain particles of the amorphous alloy Fe l _ x C x 
(x= 0.2-0.3) (ref. 57). Details regarding all tested systems are 
collected in the Supplementary Table SI. The presented analysis 
is focused on the primary relaxation time t(7), because the 
modern broad-band dielectric spectroscopy can yield sets of high 
precision and consistent experimental data 3 ' 13 . For the 
presentation of experimental data, we apply the linearized 
derivative-based analysis introduced in Drozd-Rzoska and 
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Rzoska 44 . The VFT equation allows the following transformation 
of t(T) data: 



dlnr 



d(i/n 



(D T T 0 ) 



-1/2 



T 0 (D T T 0 )- 1/2 



= B- 



For the critical-like equation (3): 

T 2 T-T c 



= AT — B 



A 
T 
(12) 



(13) 



where H'JT) stands for the normalized apparent enthalpy. 

The linear behaviour of equations (12) and (13) show the 
valid temperature domain of the VFT and 'critical-like' 
equations. The subsequent linear regression fit can yield 
optimal values of relevant parameters 44 . Results of such 
analysis for selected systems are shown in Fig. 2a-e. It is 
visible that a clear 'critical-like' parameterization is preferred for 
all LC rod-like compounds, as well as for sorbitol, ethanol and 
polystyrene (Fig. 2b). This behaviour can be linked to the 
uniaxial structure of molecules, which may lead to the local 
orientational symmetry (OS) arrangements 28 . The 'critical-like' 
behaviour also dominates the dynamics of plastic crystals 
(ODICs) 28 ' 29 , where a positional symmetry (PS) is a 
fundamental feature (Fig. 2a). Regarding EPON 828, glycerol 
and DEP, the critical-like parameterization definitively fails 
(Fig. 2c). Because these compounds cannot be linked neither to 
OS nor to PS cases, we assigned them as no-symmetry (NS) 
compounds. 

For 'NS' glass formers, the VFT equation fits the experimental 
data fairly well as shown in Fig. 2e. It is noteworthy that for 
ethanol and sorbitol, both equations (1) and (3) seem to offer 
the acceptable parameterization (Fig. 2b,d). 

For a further insight into the optimal t(T) parameterization, 



an analysis exploring I£q{T) behaviour has been carried out 
(Fig. 3 and equation (11)). It is visible that for none of the tested 
compounds, the evolution expected for BA 49 U^q(T) = const) 
or WM 19 (& = /- o 1 (r = 0)=0) relations take place. Solely, 
patterns expected for the VFT or the 'critical-like' description 
occurs. Experimental results show that such patterns are 
associated with three 'terminal, convergence' cases: (i) 
n= -(\/b)= -1/ /- o 1 (T=0) « 0.17, for ODICs and SGL 
(PS case); (ii) n«1.53, for all LCs and also low-molecular- 
weight sorbitol, propanol and polystyrene (OS case); (iii) nsil, 
for DEP, EPON 828 and glycerol (NS case). Worth to mention is 
that the points (i) and (ii) are related to approximately the same 
value of the 'critical' exponent 0 = 8.5~12 in equation (3), 
being remarkable that the I^(T) plot allows to separate OS and 
PS cases. 

A summary of the analysis covering 16 glass-forming systems is 
presented in Fig. 4 and in the Supplementary Table SI. Apart from 
the already-mentioned manifestation of three 'characteristic' values 
of the coefficient «, which can be related to the local symmetry, two 
more issues are notable. First, the clear access of the SGL system to 
the PS group, together with ODICs. Second, the 'intermediate' 
cases between three characteristic values of the coefficient n. For 
ethanol, the coefficient is nxl.2, whereas for propanol, being a 
'longer' and thus more uniaxial compound and from the same 
homologous series, the coefficient is nxl.5, as for LCs and 
polystyrene. One of the 'intermediate' cases is cyclooctanol (plastic 
crystal), for which nx0.8. In this case, the discrepancy may be 
linked to some distortion from the optimal positional ordering. 

Consequences for the configurational entropy. Assuming the 
validity of equation (4), recalling equation (18) from Martinez- 
Garcia et at and the definition of the DO index by equation (8), 
the following relation can be formulated: 



Jdo(T) 
(T/T) 



dln(£a(r)) 
d{l/T) 



1 dSc(T) 
S C (T) d(l/T) 



(14) 
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Figure 2 | Linearized derivative-based presentation of t(T) data. Results of the linearized, derivative-based transformations 44 of x(T) experimental data, 
with H' a = H a /Rand H a = rJIn i/d(1/T). (a-c) Refer to (equations (3) and (13)). (d,e) Refer to the VFT dependence (equations (1) and (12)). The domains of the 
validity of the given description are shown via the linear behaviour. The blue dashed lines in a-c stand for 'terminal critical exponents' cj> = 8 and 12, respectively. 



By applying experimental findings from the previous section, the 
above equation can be written as: 



dSc(T) 
d(l/T) 



+ 



t + b/T 



S c (T)=0 



The latter directly yields the following solution: 



Sc(T) = S 0 







- l/b 




H 


)©] 


= S 0 




T 



(15) 



(16) 



where the 'singular' temperature can be defined as T N = \ b/a\ and 
the power exponent as n = —l/b= — l/(/ I J 0 1 (r= 0)). 

Equation (16) represents a new generalized dependence for the 
evolution of the configurational entropy of glass-forming liquids, 



allowing to construct the normalized/scaling plot shown in Fig. 5. 

The scaling Fig. 5 covers the variety of ultraslowing glass- 
forming systems ranging from low-molecular liquids and 
polymers to plastic crystals, LCs and SGLs. Clearly visible are 
three characteristic (universal?) values of the power exponent: 
n x 1.53, n=l and ««0.17, which can be linked to different 
dominating molecular symmetries, as shown in the right part of 
Fig. 5. Worth to mention is also that equation (16) has been 
derived by the solution of the first-order differential 
equation (15), without underlying assumptions that were used 
in the derivation of equation (5). The latter is recovered for n= 1. 
It is also noteworthy that equation (16) indicates a possible link 
between the 'critical' temperature (7c) in equation (3) and the 
VFT singular temperature (T 0 ) in equation (1)). 
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Figure 3 | Reciprocals of the experimental DO temperature index. The temperature evolution of the inverse DO index for tested glass-forming materials 
given in the figure. The analysis disentangles three representative cases: (a) the 'PS' (ODICs and SGL); (b) the 'OS' (LC, propanol and polystyrene); 
and (c) the lack of dominant symmetry (NS). Note the convergence emerging for 1// DO (T=0). 
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Figure 4 | Coefficient n versus the exponent <j). Experimental dependencies of the symmetry coefficient (n) versus the power exponent <f> (equations (3) 
and (13)). For systems with the dominant VFT parameterization the exponent has been estimated as the terminal case 0(T->T g ). Convergence values of 
the coefficient (n) related to VFT or 'critical-like' behaviours are indicated. The separation of OS and PS cases indicates the role of the symmetry. 



Discussion 

In the year 2008, Hecksher et al} 5 boosted the discussion 
regarding one of the most fundamental experimental artifacts 
within glass-transition physics, namely the validity of the VFT 
equation, the existence of the 'final temperature divergence' < T g , 
and the ultimate way of t(T) or t](T) parameterization. This 

6 



discussion and subsequent papers • ' have indicated that 
parameterizations not associated with a 'finite temperature 
divergence' may be crucial. All these findings indirectly 
questioned the connection between the glass transition and a 
hidden phase transition below T g . However, parallel studies 
exploring the linearized derivative and distortion sensitive 
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Figure 5 | Scaling configurational entropy plot. It shows the experimental evolution of the normalized structural entropy for selected experimental 
ultraslowing systems. The analytic evolution related to PS, NS and OS patterns (solid curves) is coupled to different 'universal' values of the coefficient 
exponent n= — 1//do(7"= 0). The related symmetries are also illustrated. 



analysis led to a clear proof that the 'critical-like' description with 
a finite temperature divergence may be optimal for liquid 
crystalline glass formers, plastic crystals (not tested in Erwin 

„27,: 



and Colby ) and even for selected low-molecular liquids and 
polymers ,2> 



This paper shows that the 'finite temperature divergence' in 
ultraslowing glass-forming systems on approaching T g may be 
feasible not only for low-molecular-weight liquids and polymers 
but also for LCs, plastic crystals and even SGL systems. 
Surprisingly, a consistency seems to be possible, despite the fact 
that a single equation for portraying the dynamics in the 
ultraslowing domain for an arbitrary glass former may not exist. 
It is worth stressing that the optimal parameterization via the 
VFT equation should be considered solely for selected com- 
pounds, and other systems such as parameterization is inherently 
non-optimal. This fact ought to be taken into account when 
validating different equations for portraying the t(T) or fy(T) 
experimental behaviour. 

We would like to stress the results of the analysis focused on 
BA 48 ' 49 and WM (MYEGA) equations 1 U9 - 26 without the finite 
temperature divergence. 

Worth to mention is also a group of further equations of such 
type 2,3,13,15 , which should be tested in the future. 

In conclusion, the obtained key forms for t(T) or t](T) 
parameterization can be linked to predefined values of the 
parameter/ exponent n= — l// r J o 1 (T = 0), as well as to the 
dominating local symmetry. They are clearly associated with a 
finite temperature divergence below the glass temperature. All 
these results indicate qualitatively new links between glass 
transition and critical phenomena. 



Methods 

Analysis of experimental data. The analysis of 'dynamic' experimental data 
employing higher-order derivatives of r(T) or ?/(T) experimental data is most often 
avoided due to scatter associated with the inherent experimental error. In this 
paper, this basic issue has been overcome due to the novel implementation of 
Savitzky-Golay (SG) filtering procedure introduced in Martinez-Garcia et a/. 52 



(equation (8)). So far, it was calculated from the 'general' apparent SA equation 
via -Ea(T) = RT\n(T(T) /tq), assuming for t 0 — 10 ~ 14 s, which is most often 
considered as the 'typical' or 'universal' value, at least for ultraviscous liquids. 
Notwithstanding, there is clear evidence that the value of t 0 ranges from 
10- n -10- 16 s (ref. 49). The erroneous value t 0 for a given experimental system 
can notably influence estimations of E a (T). In this paper, experimental values of 
EJX) were calculated in a novel way, avoiding T 0 -related problems. The AAE 
was calculated via H a ( T) = H A ( T) /R = d In t/d( 1/ T) . Ultimate estimation was 
supported by the SG filtering procedure. Subsequently, the following 'general' 
first-order differential equation was solved: 



BKWf) , E 'MT) 



(l/T) 



ap/T) (i/T) 
{(l/T) -> oo: E a (0),E a (T N ) -> oo} 
{(l/T)->0:E a = cte} 



(17) 



The 'noise-related' scatter of the calculated E a (T) values was minimized by the SG 
filtering procedure for each tested glass-forming system. In a next step, the DO was 
calculated via the equation (8). 

Derivations of the DO temperature indexes for selected equations describing 
t(T) or r/(T) evolutions. The form of I DO (T) for VFT equation (1) and BA relation 
t(T) = To exp(A/T D ) can be found in Erwin and Colby 17 . The DO temperature 
index for the WM (MYEGA) equation can be obtained by linking equation (2), the 
SA equation and 7do(T) definition (equation (8)): 



JSa(T) = JUCexpl - J , andthenl^ 1 = 



Km or =-ra(c/q = ? 

(18) 



Forms of IuoiT) for the critical-like case in the ultraslowing/ultraviscous dynamic 
domain can be obtained by linking equation (3) and the general SA dependence 

T(r) = T 0 exp(4(T)/Rr) : 



E i (T) = RT 



-<t>]n 



T- T c 



hi 



= RT 



(T-T C )T 



+ In- 



to 

(19) 



E,(r)=<(r) /r=t&(~^t) + t1ii (y) +ct (20) 

This work focuses on the ultraviscous/ultraslowing dynamic domain close to T g , 
that is, for T— » T c , hence: 



Determination of the AAE from r(D experimental data. For estimating the DO 
temperature index (/do(T)), the experimental evolution of the AAE is needed 



£ a (T) sa Tin 



T-Tc 



(21) 
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Figure 6 | Analytic and approximated evolution of the 'critical-like' DO temperature index. Evolutions of /go -(analyt.) (black curve), / D r o (approx.) (blue 
curve) and their reciprocals (right side), compared with the empirical DO index for C7-OH (ODIC, plastic crystal) where the critical-like t(T) evolution is 
optimal. 



The latter equation can be re-written as: 



£ a (T) = 7T c ni 



T-Tc 



*/Tc 



: TT C S(T) 



(22) 



The function 3(T) = \n[T/(T - Tc)f' c can be expanded into the Taylor series 
because of its inherent features (<j>/Tc « 1) if \S(T) — 1| < 1. This condition is 
valid for ultraslowing ODICs, LCs and SGL systems tested in this paper. Taking 
this into account, one obtains: 



00 ( _ i y +1 

E 3 (T) = TT C J2- — 



\T-T C J 

Then, recalling 7 DO (r) definition (equation (8)): 

Wt)=i§c!t) = 



TTr 



T 



T-Tc 



4>/Tc 



( 










-In 









T-T c 



(23) 



(24) 



The evolutions of I$q(T) and lfl^(T) are presented graphically in Fig. 6. 

The analytic form of the index is given in the bracket in equation (24) 
and shown graphically in Figure 6, together with the approximation 
^do = ^DO " (approx.) = (f>/(T — Tc). This figure is supplemented by the 
experimental DO index for C7-OH (ODIC; plastic crystal) where the critical-like 
parameterization of z(T) is optimal (see Fig. 2). The validity of /^q (approx.) for 
portraying experimental data in the ultraslowing domain (below T B ) is clearly 
visible. 
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